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FOREWORD 


Significant  differences  in  underwater  shock  wave  peak  pressure  and  decay 
constant  have  been  measured  when  identical  explosive  charges  have  been  fired  at 
different  test  sites.  This  "Key  West  Effect"  has  resulted  in  uncertainties  in 
interpreting  explosion  test  data.  This  report  presents  a  method  for 
quantitatively  describing  these  and  other  changes  in  underwater  explosion  shock 
wave  parameters  in  terms  of  the  sound  speed  and  density  of  the  ambient  water. 

This  work  is  part  of  the  Explosives  Development,  Effects,  and  Safety  Program 
of  the  Naval  Sea  Systems  Command  and  was  supported  by  Task  Area  SF -33-354-391 . 


F.  PROCTOR 
By  direction 


Accession  Jor_ 

1st  IS  GRAM 

DTIC  TAB 

Unannounced 

Justification — . 


Distribution/  _ 

~  Availability  Codes_ 
~  Avail  and/or 

iDlSt  1  Special 


U1UV 


1/2 


NSWC  TR  80-491 


CONTENTS 


Section  Page 


1  INTRODUCTION .  5 

2  SIMILITUOE  EQUATIONS .  12 

Generalized  Similitude  Equations .  14 

3  RESULTS .  16 

4  DISCUSSION .  17 

APPENDIX .  A-l 


ILLUSTRATIONS 


Figure 


Page 


1.1 

1.2 

1.3 

1.4 


Generalized 
Generalized 
Generalized 
Generalized 
Density. . 


Similitude  Equation 
Similitude  Equation 
Similitude  Equation 
Similitude  Equation 


Scaling  of 
Scaling  of 
Scaling  of 
Scaling  of 


Peak  Pressure. . . . 
Decay  Constant. , . 

Impulse . 

Energy  Flux 


6 

7 

8 

9 


TABLES 


Table  Page 

1.1  Key  West  vs  Dahlgren  Differences  (Average  Values  for  2 

Sites) .  11 


3/4 


NSWC  TR  80-491 


SCALING  UNDERWATER  EXPLOSION  SHOCK  WAVES  FOR  DIFFERENCES 
IN  AMBIENT  SOUND  SPEED  AND  DENSITY 


1.  INTRODUCTION 

In  the  measurement  of  underwater  explosion  shock  waves  it  has  generally  been 
assumed  that  local  ambient  water  conditions  have  a  negligible  effect  on  the 
recorded  pressure  signature.^  The  exception  to  this  is  when  distortions  due  to 
refraction  occur.  A  recent  set  of  data  shows  that  even  with  no  refraction,  the 
shock  wave  pressure  signature  may  be  strongly  dependent  on  the  temperature  and 
salinity  of  the  water. 

In  an  underwater  explosions  field  test  conducted  by  the  Naval  Surface  Weapons 
Center  (NSWC)  the  test  site  had  to  be  moved  from  one  location  to  another  in  mid¬ 
program.  After  a  number  of  shots  had  been  fired  in  the  Potomac  River  at  Dahlgren, 
Virginia,  in  January  of  1979,  the  operation  was  suspended  due  to  environmental 
considerations.  The  remaining  shots  were  fired  in  the  Atlantic  Ocean  near  Key 
West,  Florida,  in  June  of  1979. 

The  program  consisted  of  a  number  of  different  explosives  fired  in  100-lb 
charges.  Duplicate  shots  of  four  different  charges  were  fired  at  the  two  sites. 
There  was  a  consistent  pattern  of  variation  in  measured  pressures  between  the  two 
sites.  This  is  shown  by  the  pairs  of  curves  on  the  left  hand  side  of  Figures  1.1 
through  1.4.  These  curves  are  the  familiar  power  law  fits*  made  by  the  method  of 
^Cole,  Robert  H.,  1948,  "Underwater  Explosions,"  Princeton  University  Press 

Reference  1 ,  p  238 
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FIGURE  1.1  GENERALIZED  SIMILITUDE  EQUATION  SCALING  OF  PEAK  PRESSURE 
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-  DAHLGREN 
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FIGURE  1.3  GENERALIZED  SIMILITUDE  EQUATION  SCALING  OF  IMPULSE 
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FIGURE  1.4  GENERALIZED  SIMILITUDE  EQUATION  SCALING  OF  ENERGY  FLUX  DENSITY 
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least  squares  to  the  experimental  data  for  peak  pressure,  decay  constant,  impulse 
ana  energy  flux  density  from  four  different  types  of  explosive  charges--each  fired 
at  the  two  sites.  They  show  that  in  moving  from  Dahlgren  (solid  curves)  to  Key 
West  (dashed  curves)  the  peak  pressure  increased  about  1 2 c_- ,  the  decay  constant 
decreased  about  1 4“' ,  the  impulse  decreased  about  4i;,  and  the  energy  flux  density 
was  not  significantly  changed.  Differences  between  the  Key  West  and  Dahlgren  test 
conditions  are  summarized  in  Table  1.1. 


2 

Gaspin  and  Lehto  have  shown  that  this  "Key  West  Effect"  is  real  and  that  it 

can  be  approximately  accounted  for  by  differences  between  the  two  sites  in  the 

2 

sound  speed,  c ,  and  the  bulk  modulus,  =  cQco,  where  o  is  the  ambient  water 
density.  From  an  approximate  physical  model  they  estimated  the  effect  on  peak 
overpressure,  PMAX,  and  on  decay  constant,  6,  at  a  fixed  distance  from  the  charge 


to  be 


pmax2 

PMAX^ 


(14,500  >  PMAX  >  1450  psi) 


(1.1) 


_2 

ei 


(1.2) 


Equation  1.1  is  in  good  agreement  with  the  experimental  data  shown  in  Figure  1.1, 
while  Equation  1.2  partially  corrects  for  the  change  in  decay  constant  shown  in 
Figure  1.2.  Lehto  and  Gaspin  do  not  give  results  for  the  effects  on  impulse  and 
energy  flux  density. 

The  purpose  of  this  note  is  to  demonstrate  the  possibility  of  accounting  for 
changes  in  the  peak  pressure,  decay  constant,  impulse,  and  energy  flux  densi  ty 

7 - : - 

GasDin,  J.  B.  and  Lehto,  D.  L. ,  "Effects  of  Water  Temperature  and  Salinity  on 
Underwater  Explosion  Shock  waves,"  NSWC  Report  to  be  published 


10 


NSWC  TR  80-491 


Table  1 

.  .1  Key  West  vs 

Dahlgren  Differences 

(Average  Values 

for  2  Sites) 

Oahlqren 

Key  West 

Change 

Water  Temperature  (Deg.  C) 

2.3 

27.8 

+  25.5 

Salinity 

7.6 

36.5 

+  23.9 

Sound  Speed  (ft/sec) 

4668 

5048 

+  8.1% 

Water  Density  (gm/cc) 

1.006 

1.024 

+  1.8% 

Bulk  Modulus  (lb/ft2) 

42.54  E6 

50.66  E6 

+19.1% 

Peak  Pressure 

— 

— 

+  12% 

Decay  Constant 

— 

— 

-14% 

Impulse 

— 

— 

-4% 

Energy  Flux  Density 

— 

— 

0% 
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(PMAX,  6,  I,  and  E)  due  to  variations  in  cQ  and  cQ  by  simply  rederiving  the 
dimensionless  variables  used  to  formulate  the  underwater  explosion  shocx  wave 
similitude  equations--and  not  discarding  c-  and  cQ.  The  pairs  of  curves  shown  on 
the  right-hand  side  of  Figures  1.1  through  1.4  show  the  results  of  this  effort. 

2.  SIMILITUDE  EQUATIONS 

The  similitude  equations  for  the  underwater  explosion  shock  wave1 


/,,l/3\a 

PMAX  =  k  f  ) 

(2.1) 

e  =  4  w1/3  ^ 

i4 

(2.2) 

I  ■  m  U1'3  ( 

■H 

(2.3) 

£  =  „  K1'3  ( 

W1/3\6 

(2.4) 

are  specific  expressions  (power  curve  fits,  y 

*  axb)  of  a  general 

scaling  expressed 

b“  , . 

^  *>  '  f(^75 

’  '  7^) 

(2.5) 

We  now  proceed  to  rederive  Equation  2.5  from  dimensional  considerations  but 
retaining  the  dependence  on  oQ  and  cQ.  The  significant  variables  for  the  problem 
are: 

p  «  p(t)  =  shock  wave  overpressure 
R  8  radial  distance  from  center  of  charge 
t  8  time 

Yq  =  total  chemical  energy  released  upon  detonation 

=  ambient  water  density 
o  J 

cQ  =  ambient  sound  speed. 
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Thus,  we  are  in  effect  assuming  that  the  overpressure,  p,  is  a  function  of  R,  t, 

Y  ,  o  ,  and  cQ  or  that 

f  (P.  R,  t,  Y0,  pQ,  cQ)  -  0  (2.6) 


where  f  Is  some  unknown  function.  These  six  variables  in  Equation  2.6  can  be 
expressed  in  terms  of  three  dimensions--mas$,  length,  and  time.  By  Buckingham’s 

3 

ir-theorem  ,  the  problem  can  be  reformulated  in  terms  of  three  dimensionless  groups 
or  tt ’  s  and  the  solution  must  be  of  the  form 

F  (ttj,  tt2,  tt3)  =  0  (2.7) 

A  convenient  set  of  n's  for  considering  this  problem  is  the  following* 


tr  =  P-  ■  =  B_ 
1  K 


(Vo2)1/3r  r 
*2  TT71  7173“ 
ro  '  o 


(2.8) 

(2.9) 


C0K01/3t 


(2.10) 


2 

where,  Kq  =  pQco  ,  which  Is  the  bulk  modulus  of  the  ambient  water. 

Making  use  of  this  set  of  tt's  Equation  2.7  can  be  written 

Is  G 
Ko 

where  G  is  some  unknown  function. 

^Buckingham,  E.,  1914,  "On  Physically  Similar  Systems-,  Illustrations  of  the  Use  of 
Dimensional  Equations,"  Physical  Review,  Vol.  4,  p  345 

♦The  reader  is  referred  to  References  3  and  4  for  the  method  of  determing 
convenient  sets  of  Vs. 

4 

Johnson,  Walter  C.,  1944,  "Mathematical  and  Physical  Principles  of  Engineering 
Analysis,"  McGraw-Hill,  (Chapter  IX) 


V/3r 

717 r 

,  0 


c,"3t 

0  0 


(2.11) 
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If  we  restrict  our  consideration  to  a  single  explosive,  we  can 
replace  in  Equation  2.11  by  W,  the  weight  of  explosive,  and  obtain 


p(R,t) 


'JV3') 

w1/3  / 


(2.12) 


which  is  similar  to  Equation  2.5  except  we  have  retained  the  desired  dependence  on 
and  c  .  Note,  if  we  assume  KQ  and  c0  are  constants,  then  Equation  2.12  is 
equivalent  to  the  familiar  shock  wave  scaling  expressed  by  Equation  2.5. 


Generalized  Similitude  Equations.  Note  that  the  similitude  eouations. 
Equations  2.1  through  2.4,  are  simply  power  curve  fits  to  experimental  data  using 
dimensionless  variables.  Thus,  our  generalized  similitude  equations  are  simply 


PMAX 


C  K^3  6 
0  0 


w 


T7T 


co! 

k2/3w1/3 


~77TJ7T 
Ko  w 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


The  exponents--a,  6,  v,  and  6-- are  the  same  as  in  Equations  2.1  through  2.4.  The 
coeff 1cients--k*,  a*,  m*,  and  n*--can  be  calculated  from  k,  £,  m,  and  n, 
respectively,  by  transforming  Equations  2.1  through  2.4  to  the  same  form  as  2.13 
through  2.16,  respectively.  Comparison  then  shows  that 


-  * 
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k*  =  k 

K 1  -  i 

(2.17) 

0 

i*  =  2.  C  K 

0  0 

(2.18) 

m  cfl 

s  ~{r-7T/3 

(2-19) 

n*  =  n  K0{5-2)/3  (2.20) 

Equations  2.17  through  2.20  can  be  used  to  compare  data  fits  made  using  Equations 
2.13  through  2.16  to  those  made  using  Equations  2.1  through  2.4. 


By  substitution  into  Equations  2.13  through  2.16,  we  determine  the  changes 
In  PMAX ,  9,  I,  and  E  at  constant  range  due  to  different  ambient  conditions.  They 
are  as  follows: 


(2.22) 


(2.23) 


(2.24) 
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Changes  in  the  similitude  coefficients— k,  1,  m,  n— from  one  ambient  state 
to  another  C3n  be  computed  from  Equations  2.17  through  2.20.  They  are 


A  - 


li 


(2.25) 


(2.26) 


m  ^ 

"T 


(2.27) 


6-2 


(2.26) 


3.  RESULTS 

The  pairs  of  curves  shown  on  the  right-hand  side  of  Figures  1.1  through  1.4 
have  been  calculated  from  the  corresponding  curves  shown  on  the  left  using 
Equations  2.13  through  2.20  and  the  values  for  ambient  water  density  and  sound 
speed  measured  on  the  respective  shots.  They  have  been  scaled  to  the  same  ambient 
water  condition,  pQ  =  1  gm/cc  and  cQ  =  5000  ft/sec.  For  these  tests  the  generalized 
similitude  equations.  Equations  2.13  through  2.16,  do  a  good  job  of  accounting  for 
the  observed  changes  in  all  of  the  shock  wave  parameters,  PMAX,  e,  I,  and  E.  The 
calculated  change  In  energy  flux  density  from  Dahlgren  to  Key  West  is  a  decrease  of 
about  0.1%  which  is  too  small  to  detect  experimentally  or  to  observe  on  Figure  1.4. 
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4.  DISCUSSION 


Inserting  current  NSWC  values^,  a  =  1.19,  3  =  0.26,  Y  =  0.90,  5  =  2.10,  for 
the  similitude  equation  exponents  for  cast  50/50  pentolite  into  Equations  2.21 
through  2.24,  we  get 


(4.1) 


e 

e 


2 

1 


(4.2) 


17 


(4.3) 


(4.4) 


The  result  for  peak  pressure  is  essentially  the  same  as  Lehto  and  Gaspin's  result 
(Equation  1.1).  The  result  for  decay  constant  is  apparently  more  complete  than 
their  result  (Equation  1.2).  The  results  for  impulse  and  energy  flux  density  are 
new  and  are  in  agreement  with  experimental  measurements. 


Trice,  Robert  S. 


MSUC  limited  distribution  report 
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In  the  final  analysis,  the  value  of  these  generalized  similitude  equations 
will  only  be  demonstrated  by  extensive  use.  This  should  include  incorporation 
into  the  shock  wave  data  analysis  programs  used  by  the  Naval  Surface  Weapons 
Center.  The  equations  needed  to  do  this  are  derived  in  the  Appendix. 
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APPENDIX  A 

POWER  LAW  FITS  TO  MEASURED  DATA 


Since  self-consistent  units  of  measurement  are  not  used  in  underwater 
explosion  research  in  this  country,  it  is  not  recomnended  to  make  power  law  fits 
to  experimental  measurements  using  Equations  2.13  through  2.16.  Instead,  we  can 
make  use  of  the  dimensionless  variables  implied  in  Equations  2.13  through  2.16 
and  scale  our  measured  experimental  quantities--R,  PMAX,  0,  I,  E--to  some  standard 
condition,  e.g.,  fresh  water  at  some  standard  temperature  and  pressure*.  Thus, 


using  Equations  2.13  through  2.16 

R* 


C/3* 

,1/3 


;1'3r 

s _ s 

wl/3 


(A.l) 


(A. 2) 

(A. 3) 


(A. 4) 


*The  author  recommends  that  a  standard  condition  be  adopted  for  reporting  measured 
values  of  similitude  coefficients  and  suggests  that  fresh  water  at  15°C  and  1 
atmosphere  pressure  would  be  suitable.  With  this  choice  p  =  0.99913  gm/cc  and 
c$  =  4809.3  ft/sec.6  s 


Sel  Grosso,  V.  A.,  "New  Equation  for  the  Speed  of  Sound  In  Natural  Water;,  (With 
Comparisons  to  Other  Equations),"  J.  Acoust.  Soc.  Am.,  October  1974 


A-l 
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where  denotes  a  dimensionless  variable  and  subscript  "s"  denotes  standard 
condi tion. 


Solving  Equations  A.l  through  A. 5  for  R$,  PMAXg,  e  ,  I  ,  Eg,  respectively, 


we  get 


■ft) 


1/3 


(A. 6) 


PMAX,  =  ~  PMAX 
S  Ko 


(A. 7) 


(A. 8) 


(A. 9) 


(A. 10) 


The  derived  quantities  Rs>  PMAXs,  @s,  I$  and  Eg  can  be  used  to  make  the  power 
law  fits  indicated  by  the  ordinary  similitude  equations  (Equations  2.1  through 
2.4);  and  thereby  obtain  coefficients--k  ,  St  ,  m  ,  n  — and  exponents— a  6  ,  y  , 

J  O  p  J  O  W  -> 

6s--scaled  to  the  desired  standard  condition.  This  can  be  verified  by  substituting 
Equations  A. 6  through  A. 10  into  Equations  2.13  through  2.16  to  get 


PMAX$ 


Kg  T 


(A. 11) 


.  J*  /  Rs  \  B 

?75  -  [rpi) 


(A. 12) 


A-2 
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I$  m*^2'^73  I  R$ 
W^3  cs  \W1/3/ 


(A. 13) 


(A. 14) 


Substituting  Equations  2.17  through  2.20  into  Equations  A. 11  through  A. 14, 
respectively,  we  get 


(A. 15) 


(A. 16) 


(A. 17) 


(A. 18) 


which  is  our  desired  result. 
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